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Abstract 

The pressure- velocity correlation and return to isotropy term in the Reynolds stress 
transport equation are analyzed using the Yakhot-Orszag renormalization group. The per- 
turbation series for the relevant correlations, evaluated to lowest order in the e-expansion 
of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean ve- 
locity gradient and its transpose. Formal lowest order Pade approximations to the sums of 
these series produce a fast pressure strain model of the form proposed by Launder, Reece, 
and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, 
the model constants are computed theoretically. The predicted Reynolds stress ratios in 
simple shear flows are evaluated and compared with experimental data. The possibility 
is discussed of deriving higher order nonlinear models by approximating the sums more 
accurately. 

*Work funded under NASA Cooperative Agreement NCC3-233. ICOMP Program Director, 
Louis A. Povinelli. 



Introduction 


The modelling of the pressure correlation and return to isotropy term in the Reynolds 
stress transport equation remains controversial. 1,2,3 Models will be derived here using the 
Yakhot-Orszag renormalization group 4 partially along the lines of our previous work 5 . The 
result is a model for the fast pressure-strain term of the form proposed by Launder, Reece 
and Rodi 6 (LRR) and a model for return to isotropy of the form proposed by Rotta 7 with 
theoretically computed constants in good agreement with accepted values. As usual in 
investigations of this sort, the priority of Yoshizawa in deriving a pressure strain model 
analytically 8 must be noted. 

The analysis requires some new ideas in renormalization group theory recently intro- 
duced by Yakhot et al 9 . As Yakhot et al 9 emphasize, the application of the renormalization 
group mode elimination formalism to shear flow creates a double perturbation series in 
powers of e, the parameter of the isotropic theory, and in powers of a dimensionless strain 
rate, r) = SK/e, where K denotes the turbulence kinetic energy, e denotes the dissipation 
rate, and S' is a measure of the mean strain: in Ref. 9, S 2 — The 

present analysis also leads to double expansions of this type, with the powers S n replaced 
by tensors homogeneous of degree n in the mean velocity gradient VU and its trans- 
pose. It will be convenient to retain the terminology of Ref. 9 and call this expansion the 
jf-expansion; when the distinction is pertinent, the expansion of Ref. 9 will be called a 
scalar ^-expansion. 

The heuristic program of evaluating all scalar amplitudes to lowest order in e has 
proven successful in the past: apparently, the e-expansion is an asymptotic series with 
sum given very nearly by its first term 10 . Unfortunately, there is no analogous basis for 
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truncating the 77-expansion. There are fundamental reasons for this distinction between 
these expansions. The present 77-expansion is tensorial: successively higher order terms 
do not introduce merely numerical corrections, but increasingly complex asymmetries into 
the theory. Truncation therefore imposes a possibly inappropriate symmetry or other 
constraint on the model. Thus, in Ref. 5 the 77-expansion for the Reynolds stress r was 
truncated at second order as suggested by previous work of Yoshizawa 11 and Speziale 12 . 
Although this type of modelling permits unequal normal stresses in a simple shear flow, 
it is not maximally asymmetric: for example, in a flow with mean velocity components 
Ui(x i,X2), a cubic model including a term t ~ VU 2 VU T +'VUVU T2 would permit nonzero 
T2 3 in the presence of vanishing dUildx%, an effect which cannot be ruled out in advance. 

Although generalizations 13 of the Cayley-Hamilton Theorem limit the number of in- 
dependent tensors S^, anisotropy and asymmetry cannot exist at all without some terms 
of higher order in 77; indeed, truncation at lowest order in 77 just produces a theory of 
isotropic turbulence. But the series truncated at any higher order can be unsatisfactory in 
flow regions in which some components are large. In such regions, the truncated 

series is dominated by its highest order terms. For the quadratic stress models of Refs. 5 , 
11, 12, this domination can produce negative normal stresses in the buffer layers of wall 
bounded flows. Increasing the order of truncation obviously exacerbates this problem. 

It follows that finite truncation of the 77-expansion is theoretically unsatisfactory. 
Yakhot et al 9 therefore propose that this expansion must be summed, even if only ap- 
proximately, and have suggested a prototype summation in a different context. It should 
be noted that the same issues arise naturally in Yoshizawa’s formalism, which also gen- 
erates infinite series in the mean velocity gradients (and in other quantities as well) for 
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correlations of interest in turbulence modeling. Yoshizawa has concluded independently 
that summation of this series is essential and has also derived a Reynolds stress transport 
model by introducing such summations 14 . 

In this paper, the perturbation series which the Yakhot-Orszag renormalization group 
generates for the correlation 


rr , dp dp 

/7u = - ( “ i aJ- +u '^ 


(i) 


is summed by a low order Pade approximation. Coefficients axe evaluated to lowest order 
in the e expansion, but the summation includes effects of all orders in 77. The result is 
essentially identical to the “model 1 ” proposed by Launder, Reece, and Rodi 6 . An entirely 
analogous treatment of return to isotropy yields a model of the form proposed by Rotta 7 . 
Combining these models leads to a preliminary Reynolds stress transport model. The 
problem of closing the Reynolds stress diffusion terms is addressed. This problem also 
leads to an infinite sum. 

While it is encouraging that renormalization group methods can be used to derive 
familiar models, the goal of this investigation is not limited to providing theoretical jus- 
tification for the LRR and Rotta models, which although widely applied are nevertheless 
deficient in several well-documented respects 1,2 ’ 3 . Instead, renormalization group methods 
together with approximate summation of the 77-expansion can be used to derive higher 
order and nonlinear corrections to these models in a systematic fashion. Explicit develop- 
ment of such models is left to future investigations. 


I. Analysis of the Pressure Correlation 
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The analysis will follow Yakhot and Orszag’s derivation of turbulence transport models 


by renormalization group methods 4 . The equation for velocity products is 


d 


duj 


du{ 


dt UiU i + u r [ Ui dx p +Uj dx p 


= -( u 'l; +u '^) +,/oV 2 ui “' _ 2 ‘' 0 i;^ (2) 


where u 0 denotes the kinematic viscosity. The product —(uidp/dxj+ujdp/ dx{) on the right 
side of Eq. (2) will become the correlation Ilij defined by Eq. (1) following elimination of 
all fluctuating modes. Its Fourier transform is 


J Uj(k - q)iq t q 2 (q p - Q p )u q (q - Q)Q q u p (Q)dq dQ /(2ir) 2d+2 
+ J Ui(k - q)iqj q~ 2 (q p - Q P )u q (q - Q)Q q u p (Q)dq dQ/( 27r) 2d+2 (3) 

Here the standard notation 

k = (u>,k) fc 2 = k-k 

is used and d = 3 is the number of space dimensions. Introduce an ultraviolet cutoff 
Aj, of the order of the inverse Kolmogorov scale (t'o/s) 1 ^ 4 ; only inertial range scales with 
k < Ad will be treated explicitly. Introducing a parameter r initially near zero, partition 
wavenumber space 0 < A: < A d into the two intervals 0 < k < A^e r and A^e < k < A^- 
Denote velocity components with wavevectors in the first interval by the superscript < 
and those in the second by >. Introducing this decomposition into Eq. (3) produces eight 
terms; however, as in analogous calculations in Ref. 4, only three will contribute at the 
lowest order in e: 
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1 = J u>(k - q)i,,q-\, p - Q p )u>(q - Q)Q,u}(Q)dqdQ/(2wf d+2 
II = J - q)iq iq ~ 2 (q r - Q r )uf(q - Q)Q,u>(Q)dqdQ/(2n) 21+2 
III = J uf(k - q)iq iq -\q p - Q r )u>(q - Q)Q,u>{Q)dqdQ 

The > modes are to be eliminated from these expressions by iterated use of the 
randomly forced Navier Stokes equations 


(-* w + I/O k 2 ) Ui (k) = -^ A 0 P imn (k) J u m (k — q ) u m ( q ) dq/( 27i) d+1 + f, 

where 

-P.mn(k) = k m Pj n (k) + k n P im ( k) 

P in(k) = 8in - ki k n /k 2 

and the Gaussian random force /, is defined by its correlation function 

(fi(k) fj(k')) = 2 (2ir) d+1 D 0 8(u + J ')S(k + k' ) 

The choice y — d generates a Kolmogorov inertial range. A detailed exposition of this 
procedure can be found in Ref. 4. It will suffice to note here that the result of the 
mode elimination is a series in powers of u < and Ao- At each order in u < , perturbation 
theory will produce a finite number of types of terms with amplitudes given as series in 
Ao- Under iterated mode elimination to the limit r — > oo, the expansion in Ao proves to 
be an expansion in powers of e = 4 y — d. This is the e-expansion of the Yakhot-Orszag 
theory. Previous experience 4 and preliminary analysis 10 suggest that the amplitudes are 
best evaluated at lowest order in e with e set to zero. This procedure will be followed here. 
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Thus, the perturbation series will be written as 

n = To + Tx + • ■ • 

where T n is of order n in u < and all amplitudes are evaluated to lowest order in e. To 
lowest order in e and SK/e 

III = j 6(q)qi q p q ~ 2 Q q Uj(k)u q (-Q)u p (Q) dq/(2n) d+1 = 0 

since incompressibility implies Q q u q = 0. Now term III is formally proportional to qi uj 
and since only indices i and j are uncontracted, the combination Q p u p or Q q u q must 
occur at all orders. Accordingly, term III vanishes identically. 

Lowest order analysis of term I gives 


j Hip q-*- v { |G(9 )| 2 Pj, (q) 2 D. dg/( 2 ^) J+1 } • ik, n r (k) 

Adding the corresponding term from II and the result of ij index interchange leads to 


Tx = T i ( 


dUf 

V dxj 


+ 


m \ 

dxi ) 


where Ti satisfies the recursion relation 


dTi _ 6 V 
"dr" “ 15 vA 2 


In Eq. (5), y has been set equal to d to obtain Kolmogorov scaling and 


( 4 ) 


( 5 ) 


V = 2D 0 S d /{2^) d 
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A = A<*e~ r 


where Sd denotes the area of the d - dimensional sphere. Integrating the recursion relation 
(5) in the high Reynolds number asymptotic limit 4 r — > co, 



( 6 ) 


Eqs. (4) and (6) give 



in agreement with the analysis of Crow. 16 
At the next order in SK/e, 


I = f Wp r 2 |iW-q) <?(-$) m-Q )\ 2 X 

Pqs (Q - q) |Q - q|- y 2D 0 • v r (k - Q)^dq/(2n) d + l ■ iQ q u p (Q ) + (ij) (8) 

where (ij) indicates the result of ij index interchange in the previous term. Adding the 
contribution from II, 


_t 1 \^ ■ du f\ du f , o (^f ■ du<_\du^}^ 


105 


K UU-' OU„\0U> /UUj OU^\OU^lW) 


where (0) indicates the deviatoric part: it is immediate that the ij index contraction of 
Eq. (8) vanishes. The amplitude T 2 satisfies 


dt 2 _ V 

dr 4i/ 2 A 4 


In the high Reynolds number limit, 
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1 1 40 3 

T 2 — — — — v 

2 105 4 3 4 


r /dUi 

dU p ' 

dUj A , 

( dUi 

dU p ^ 

dU p ] 

1 

00 

N2 

'qT 

H 

^9 

+ ~^ JL 
Ox, j 


ydxp 

dxi j 

dx]\ 


( 0 ) 


+ (u) 


1 

~ 

(dU, 

dl V 

dU, n , 

fdUi 

dU p \ 

dUA 

21 V 

16 | 

\dx p 

+ - 
dxi j 

lax; +2 l 

V @ x p 

+ IT*- 

OXi ) 

CD 


T CO) 


( 9 ) 


The next order will produce a term T 3 containing cubic products of velocities u < . 
In view of the form of the LRR model, it is reasonable to ask whether a term with only 
one gradient, proportional in the high Reynolds number limit to rV17 might occur at this 
order. Such terms do occur, but they cancel. Evaluation of T 3 proves to require expansions 
of the projection operators to second order, leading instead to terms S ^ homogeneous of 
degree three in the mean velocity gradient and its transpose. In general, the term T n of 
order n has the form S^ n \K/e) n . As noted in the Introduction, it will be imperative to 
include effects of all orders in SK/e in the model, but because the terms T n involve ever 
higher order derivatives of the transverse projection operators, they do not have an obvious 
law of formation. Therefore, an exact summation does not appear feasible. 

It can be verified that the terms in braces in T\ and T 2 are the lowest order terms in 
the expansion of the correlation spectrum tensor E mn with the property 




J E mn (q) dq 


and that this identification holds to all orders. Thus, perturbation theory gives the stan- 
dard result 


n „ = J, 


9 [Qip^jq + QjpEiq] dq 


dUp 

dx q 
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but with a series for the right hand side which can be evaluated explicitly to any finite 
order. In fact, in view of Eqs. (7) and (9), 



A simple approximate summation is obtained by introducing into Eq. (10) the perturbation 
series 5 for u t Uj^ in the form 


v 


dUj dUA 

dxj dxi ) 


-W” + £ s<”>(jr/ E )" 

n> 2 


and dropping the quadratic terms. The resulting model, 




_( 0 ) dUp^ 


U{Up 


( 0 ) 


dUj 


<o) dUp 


(o) dU » 

dT 

l(0) 


( 0 ) 


p J 




(ii) 


with 


Cri = ^ = - 7619 C f2 = = .0952 (12) 

agrees with the perturbation series (10) to terms of order S^. However, unlike the explicit 
quadratic model which results from simply dropping the O(S^) terms in Eq. (10), this 
model includes effects of all order in $>K/e. The consequences of this fact will be discussed 
later. This type of summation has also been applied by Yoshizawa 14 . Eqs. (11) and (12) 
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can be compared with Launder, Reece and Rodi’s “model 1”, Eq. (11) with the empirically 
adjusted constants 


Cr i= 7636 Cr 2 = .1091 (13) 

The constants Cfi and C * 2 are not chosen independently; instead, to insure certain con- 
sistency properties 6 ’ 7 , they are linear functions of a parameter C 2 . The LRR model cor- 
responds to the choice C 2 = .4; Eqs. (11) and (12) correspond instead to the choice 
C 2 = 8/21 ~ .36. 

The approximate summation used to derive Eq. (11) can be systematically gener- 
alized to generate an infinite number of models for II l] . For example, suppose that the 
perturbation series for r is introduced into the cubic terms in the perturbation series (10) 
instead of in the quadratic terms as above. This substitution will produce a model which 
can be written symbolically in the form 

n ~ S (1) + S< 2) + r(S (1) ' + S (2) ') 

where r denotes a sum of matrix products in all possible orders of r and terms S^. 
The requirement that the original series agree to order with the approximation when 
r is replaced by its perturbation series determines this approximation uniquely. 

There is obviously a strong formal resemblance between this approximation scheme 
and Pade approximation. They differ in the introduction of r as an auxiliary quantity, but 
more fundamentally in the non-commutativity of all of the variables r, Vf7, and X’U 1 . 
It will be evident in Sect. IV that introduction of these approximations into a Reynolds 
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stress transport equation will lead under suitable hypotheses to algebraic models for r of 
the form 


r ~ K (<S (1) + . . . + S (n) ) + r (<S (1) ' + . . . + S (n) ') 

which formally express r as a ratio of polynomials of degree n in V U and V U 1 . It is 
noteworthy that the LRR model is the lowest order member of this series. Moreover, the 
analysis suggests that there is no unique optimal form for the closure of 77 in terms of 
r, VU, and VU T ] instead, there is a series of approximations of (presumably) increasing 
accuracy. 

The LRR model has been criticized by Shih and Lumley 2 because it fails to insure 
realizability at the limit of two component flow. Speziale 1 has found different limitations 
when this model is applied to homogeneous shear flow. It is not clear whether a different 
summation procedure would lead to models which could answer these criticisms, perhaps 
of the forms proposed in Refs. 1 and 2. However, these references also indicate that 
what conditions a good model should satisfy is itself a somewhat controversial question. 
Accordingly, the agreement of this theory with a plausible and often used pressure strain 
model is encouraging. 

II. The Return to Isotropy Model 

The analytical description of return to isotropy is no less controversial than the mod- 
eling of the fast pressure strain term 3 . In the usual approach to turbulence modeling, in 
which correlations generated by Reynolds averaging are closed phenomenologically, this 
process is considered to result partly from the pressure correlation through a “slow” term 
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independent of the mean flow, and partly from the deviatoric part of the dissipative cor- 
relation From this viewpoint, the analysis in Sect. I is incomplete because 

it discloses only a term proportional to the mean velocity gradient, but no slow term. 
The return to isotropy will be derived here by renormalization group methods following a 
suggestion of Yakhot 17 . 

From the renormalization group viewpoint, it is natural to investigate the return to 
isotropy, even independently of the stress transport equation, by writing the perturbation 
series for 


Ui ^t + Ui 15t = / U '^ ~ + (*i) ( 14 ) 

This perturbation series differs from the perturbation series for the Reynolds stresses previ- 
ously reported 5 only in the occurence of an additional factor —ioj in all frequency integrals. 
It is therefore natural to identify the sum of this series as a functional of the Reynolds 
stresses which, like the slow term of turbulence modeling, is independent of the mean flow. 
By substituting the Navier-Stokes equations for the time derivatives in Eq. (14), one just 

recovers the equations of motion for velocity products, Eq. (2); thus, the quantity 

/ dui dui\ 

n ^ = \ Ui ~dr +Uj ~dt) 

which results from eliminating all fluctuating modes from Eq. (14) contains contributions 
from the pressure correlation through terms containing the transverse projection operator 
and contributions from the dissipative correlation through the terms containing Vq. 

By evaluating these terms in the form of Eq. (14) and insisting that the sum be 
independent of the mean flow, we are extracting properties which were not disclosed in the 
analysis of Sect. I. 
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The analysis is straightforward. Only the deviatoric terms require attention because 
the part of the correlation proporational to 6{j contributes to the transport equation for 
K which has been analyzed by Yakhot and Smith 15 . The lowest order deviator appears at 
first order in r/; to lowest order in e 


*-*(g+g) 


(15) 


where 


= (m 
dr 15 uA 2 V ; 

In view of the form of the Rotta model, it is reasonable to seek terms at the next order pro- 
portional to U{Uj. As in Sect. I, such terms do appear, but cancel exactly. This apparently 
ubiquitous cancellation was also obtained by Smith and Reynolds 18 in an analysis of the e 
transport equation. Accordingly, the second order analysis in r/ produces quadratic terms 
in the velocity gradients. Finite truncation of this series violates the requirement that 
return to isotropy be independent of the mean flow. Therefore, we must seek a reasonable 
approximate summation. The form of the lowest order term given in Eqs. (15) and (16) 
suggests 


n ‘i = I"' 


tm 3Uj 

V dxj dx 


f). 


Tx 

V 


UiU 


-(0) 

J 


Despite its triviality, this replacement does produce an approximate sum which agrees 
exactly with perturbation theory to lowest order. It therefore can be considered a type of 
Pade approximation. This lowest order summation yields 
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(17) 


= -^u7^i7 <0, 


where 


dZ _ 1 V 
dr 15 v 2 A 2 

At the infinite Reynolds number asymptotic limit, Eqs. (17) and (18) iterate to 4 


(18) 


n[j = -c R ^ur*f» (is) 

where, in the Yakhot-Orszag theory, Cr = P/e ~ 1.6. Equation (19) is therefore simply 
the standard Rotta model with Rotta constant ~ 1.6 in agreement with an earlier proposal 
of Yakhot 17 . 

A preliminary discussion of higher order summation may be appropriate. By analyz- 
ing the spectral dynamics of the return to isotropy, Weinstock 3 concluded that the shear 
and normal stresses relax at different rates. Although this behavior is obviously not ac- 
commodated by the Rotta model, it is consistent with the present theory: the perturbation 
series for II' is obtained from the series for r by multiplying the term of order n by the 
factor C n e / K for some constant C n . The C n are all unequal; therefore, the Rotta model 
is not exact. Now comparison with the series for r shows 5 that relaxation of the shear 
stress is governed by the linear term S^, whereas relaxation of the normal stresses is 
governed by the quadratic term S^ 2 ^. Since C 2 ^ Ci, these stresses relax at different rates. 
The difference is suppressed in the Rotta model, which arose in the present formalism by 
replacing all of the C n by C \ . 
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III. Reynolds Stress Transport Models 

The renormalization group describes the effect of the universal small scales of tur- 
bulence on the large scales. Therefore, convection and production, which are determined 
entirely by the large scales, cannot be derived by the renormalization group and must be 
introduced instead by Reynolds averaging 15 . Combining the usual convection and produc- 
tion terms with the pressure correlation and return to isotropy models, Eqs. (11) and (19), 
gives the Reynolds stress transport equation 


duiuf ' °) ^ y dujuj^ 


dt 

+ 


dx. 


l K ( ' d l± 

5 \dxj 


+ 


dUj 

dxi 


( dUj dU t \ W 

\ UiUp dx p +UjUp dx p ) 

+ Cr l S- + UjU^ 0) 


- Cr u,u/ 0) 

dUj\ {0) 

dx p ) 


+c «(^ <0, ^ + ^ <0, © <0)+diff - (20) 
where “diff” denotes diffusion terms which will be discussed later. The predictions of 
models of the form (20) for homogeneous shear flow, in which the diffusion terms vanish, 
have been analyzed definitively by Speziale. 19 It will be useful to generalize this analysis 
somewhat and consider any simple shear flow with exactly one nonvanishing mean velocity 
gradient component S = dU\ {dxi , in which diffusion of all Reynolds stress components is 
negligible, and in which uJuj/K is constant. These are the conditions under which Rodi’s 
algebraic models 20 can be derived and include homogeneous shear flow as a special case. 
It follows from Eq. (20) that under these conditions, 


(u lU2 /K) 2 =^~ 


P]± 

15 Cr + P/s 


~!( 


Pf* 

C R + P/e - 1 


(cvi-ir 
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— 4 (Cri — 1) Cf2 + Cj 2 


l K = 
-mrf 0) /K = 


P/£ 

Cr + P/c- 1 

P/e 

Cr + P/s-1 


“g (Cfl “ 1) + g C„ 

-| (C« - 1) + | C« 


where P = — is production. Note in particular that the ratio 


( 21 ) 


— UiUi^ /u 2 u 2 ^ — 


— f(CVl — 1) + |Cf2 
-§(CVi-l) + fCV 2 


( 22 ) 


is independent of both P/e and of the model constant Cr . 

These results can be applied to model calibration. Rotta noted 7 that certain consis- 
tency conditions require 


8 Cri -Cr2=6 (23) 

In view of Eqs. (22) and (23), the model proposed here, which is defined by the values of 
Cri and Cf 2 of Eq. (12) is the unique model of its form for which 


—uiu\^/u 2 u 2 ^ = 4/3 - (24) 

The validity of Eq. (24) can be assessed from the experimental data summarized in Table 
II. Other Reynolds stress ratios can also be readily evaluated. Substituting Eq. (23) into 
Eqs. (21) 

uruT {0) /I< =c* (4 Cri - 

/ 22 
-W&V/K =C* f 10C H - y 

-u^ 0) /K =C* - y 
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(25) 


(uiw/if) 2 =~C- - j(C') 2 (33C?, -42 C n + 13) 


where 


C* = p / £ 

C* + P/e- 1 

Simple shear flow data suggest that ujup 0 '/!^ —v^ui^/K, and —u^uz^/K are all 
positive. In view of Eq. (25), this forces CV 1 to be in the narrow interval 


— < C n < 
15 “ r 


7 

9 


from about .73 to .78. Values of the Reynolds stress ratios for various values of Cn 
appear in Table I. No value of Cn gives entirely satisfactory values for all of the ratios; in 
particular, the ratio u 1 u 1 ^°’ 1 / K ~ .4 in homogeneous shear flow cannot be obtained from 
a model of this form. 

Phenomenological modelling has suggested numerous modifications of the original 
LRRform. The simplest modification retains the form of Eq. (11), but drops the constraint 
expressed by Eq. (23). Thus, the constants Cn and Cn are considered independent, as in 
Launder, Reece and Rodi’s “model 2” 6 , or in the model of Speziale, Sarkar, and Gatski 1 . 
Justification of this step requires that the idea of modelling each individual correlation in 
the exact Reynolds stress transport equation be abandoned; instead, the model is proposed 
for the entire equation at once. A reasonable model of this type is Eq. (11) with the 
constants 


Cn = 


43 

63 


Cr 2 = — 
63 
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chosen so that the Reynolds stress transport equation reduces to the nonlinear model of 
Ref. 5 when convection and diffusion of Reynolds stress are negligible. The values of the 
Reynolds stress ratios for this model in Table I are in excellent agreement with shear flow 
data, but the fundamental validity of such modelling is unclear. An alternative procedure 
is to explore the higher order summations described in Sect. I. 

A point emphasized by both Speziale 19 and Reynolds 21 can be noted here, that while 
the elementary eddy viscosity formula and its nonlinear generalization 5 11 1,1 2 applied to 
simple shear flow give 


u^/K-SK/e u^rf 0) /K ~(SK/e) 2 (26) 

for any value of SK/e, models of the form of Eq. (20) give (26) for moderate SI\/e, but 

U 1 U 2 /K ~ const. mm^/K ~ const. (27) 

in the rapid distortion limit SK/e — ► 00 . The finite limit expressed by Eq. (27) also 
holds for algebraic Reynolds stress models derived following Rodi’s original suggestion 20 , 
for example for the model proposed in Ref. 22, 


2 K 2 /e 
UT ~ 15 C R + (P/e - 1) 

The behavior (26) is certainly incorrect when SK/e — ► 00 ; this reflects the derivation, for 
example of the relation vt ~ K 2 /e from Kolmogorov scaling. This derivation assumes 
quasi-static spectral evolution in which the Kolmogorov spectrum instantaneously adjusts 
to local conditions in both space and time. The finite limit (27) at large SK/e can only 
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occur if effects of all order in SKjz are included. No finite truncation of the 77 -expansion 
will have this behavior. 


IV. Algebraic Reynolds Stress Models 

The approximation, due to Rodi 20 , of the Reynolds stress transport equation by an 
algebraic model under the conditions of semi-homogeneous flow (negligible diffusion of r 
and rfK approximately constant) takes the form 


P-e 

K 


-(o) _ 


UiUj — ( U{Up 


dUj , dUi\ 

dr r +u ’ u ’er P ) 


(0) 


+ n'i. + n ij 


(29) 


where II and II 1 depend on r and VC/. Explicit solutions for r can be obtained, at least 
in principle, for any such approximation 23 . Briefly, one introduces a basis for polynomials 
in VC7, and VU T . The basis contains 11 terms of homogeneity order n < 5. Writing r as 
a sum of these terms with unknown coefficients and substituting in Eq. (29) leads to the 
explicit expression 


t/K = ]T H { t m) sj n) (VU,VU T ) 


(30) 


where is a scalar function of VC/ and VC/ T such that 


ffj m) ~ |VC/| m 

when |VC/| —* 00 . The assumptions made on the approximate summations require m + n = 
0; thus, t/K is bounded when SK/e — * 00 . For example, the familiar eddy viscosity 
formula is replaced in Eq. (30) by a term '*• 
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This is the type of eddy viscosity modification sought by Horiuti 24 , but with the bound- 
edness property expressed by Eq. (27). 

r 

Pope observed 23 that the coefficients H^~ n ^ in Eq. (30) would certainly be intractably 
complex; although they could be explicitly exhibited by symbolic computation, the result 
would only pertain to the particular implicit equation for the Reynolds stresses assumed 
initially in Eq. (29). Therefore, it is equally reasonable just to postulate simple forms for 
the functions for example, Eq. (28) suggests the possibility 

2 1 

15 Cr - 1 + y/UlSKIe. 

based on the identification t/K ~ This type of modeling could be particularly 

interesting when applied to the coefficients of the quadratically nonlinear models of Refs. 
5, 11, and 12. 

Some modifications of the summation procedures used here suggest themselves. First, 
it is perhaps closer to the spirit of Pade approximation to substitute 77 itself into its pertur- 
bation series instead of r. The procedure would lead to fast pressure strain models which 
are explicit function of VU and VU T . Although previous experience strongly suggests 
the appearance of r in the model, this unconventional procedure may deserve further con- 
sideration. A second related possibility is to apply the Pade method to the perturbation 
series 5 for r directly. This also will produce a family of implicit models linear in r and of 
all orders in VU and VU T , and to explicit models of the form (30). The lowest order such 
model, analogous to the LRR model, would have the form 


77 (- x > = 


u.u; = g KSi 


(dUi 
13 V \dxj 


+ 


dUj 

dx{ 


+ C\ 


K 




<o>5i +irir <o)tti 
p - + t W 


3x, 


^ K 

+ C 2 ( UiUp 

£ 1 


( 0 ) 


dl L 


dxj 


v , (o) 

- + w’a- 


) 


with model constants C\ and C?. This is actually the procedure followed by Yoshizawa 14 , 
although Yoshizawa’s two-scale perturbation theory leads naturally to a transport model 
rather than to an algebraic model. 

V. The Diffusion Term 


The diffusion of Reynolds stress arises from the triple correlation 


d l \ 

dT, {u< 

Lowest order analysis of this term will lead to a diffusive term 

d duiuj 

dx p ° r V dx p 

where a r ~ 1.4 in the high Reynolds number limit 4 , otherwise expressed, to an isotropic 
diffusivity for Reynolds stress k — a T v. 

In Ref. 22, we analyzed the diffusion of a passive scalar to second order in the 77- 
expansion, and found, as in a similar analysis by Yoshizawa 25 , corrections leading to 

( 31 ) 


, ^ - K dUi 

K{j — K €{j -f- 7C 

e dxj 
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Unlike the isotropic diffusivity, this model permits a nonzero diffusivity k \ 2 - However, it 
rules out inequality of the normal diffusivities. At the next order, perturbation theory will 
correct Eq. ( 31 ) with terms quadratic in the velocity gradients. In this theory, unequal 
normal diffusivities are possible. But now the discussion in the Introduction applies again: 
it is necessary to sum this 77-expansion. The type of Pade approximation used in Sects. 
I-II may lead to a diffusivity dependent on the Reynolds stresses, as in the passive scalar 
models of Rogers et al. 26 This type of diffusivity has also been proposed by Launder, Reece 
and Rodi 6 for the Reynolds stresses. The details are considerably more elaborate than for 
the pressure correlation, and this possibility will be left for future investigation. 

VI. Conclusions 

The present analysis of the Reynolds stress transport equation, based on the Yakhot- 
Orszag renormalization group and (tensorial) 77— expansion summation as suggested by 
Yakhot et al. 9 , has led to a model transport equation incorporating the well-known LRR 
and Rotta models. The analysis gives theoretical support both to these models and to the 
constants sometimes used with them. More significantly, it exhibits the LRR and Rotta 
models as lowest order approximations, and therefore also supports their replacement with 
higher order nonlinear models which would be deduced by more accurate approximate 
summations. The consistency of the analysis with higher order effects like the unequal 
relaxation rates of shear and normal stresses has been discussed. 
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